Gluon radiation off hard quarks in a nuclear environment: opacity expansion 



Urs Achim Wiedemann 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 
(February 1, 2008) 



O ■ 

o ■ 
o : 

<: 

in : 



(N ' 

> : 

ON . 
(N ■ 

in : 

o . 
o ■ 
o ; 

^ : 

a: 
i 



X 



We study the relation between the Baier-Dokshitzer- 
Mueller-Peigne-Schiff (BDMPS) and Zakharov formalisms for 
medium-induced gluon radiation off hard quarks, and the ra- 
diation off very few scattering centers. Based on the non- 
abelian Furry approximation for the motion of hard partons 
in a spatially extended colour field, we derive a compact di- 
agrammatic and explicitly colour trivial expression for the 
AT-th order term of the k^-differential gluon radiation cross 
section in an expansion in the opacity of the medium. Re- 
summing this quantity to all orders in opacity, we obtain Za- 
kharov's path-integral expression (supplemented with a reg- 
ularization prescription). This provides a new proof of the 
equivalence of the BDMPS and Zakharov formalisms which 
extends previous arguments to the kx-diffcrential cross sec- 
tion. We give explicit analytical results up to third order in 
opacity for both the gluon radiation cross section of free in- 
coming and of in-medium produced quarks. The AT-th order 
term in the opacity expansion of the radiation cross section 
is found to be a convolution of the radiation associated to AT- 
fold rescattering and a readjustment of the probabilities that 
rescattering occurs with less than AT scattering centers. Both 
informations can be disentangled by factorizing out of the ra- 
diation cross section a term which depends only on the mean 
free path of the projectile. This allows to infer analytical ex- 
pressions for the totally coherent and totally incoherent limits 
of the radiation cross section to arbitrary orders in opacity. 



I. INTRODUCTION 

Hard partons, produced in relativistic heavy ion colli- 
sions at RHIC and LHC, will undergo multiple rescat- 
tering inside the nuclear environment before entering 
the hadronization process outside the nuclear environ- 
ment (or in a very dilute one). This follows from stan- 
dard formation time arguments. Prior to hadronization, 
medium-induced radiative energy loss is expected to be 
the main medium modification encountered by hard par- 
tons. The size of this effect depends on the density and 
nature of the medium and may be significantly enhanced 
in a deconfined partonic medium [ 0, 0] ■ A correspond- 
ing strong medium-modification of the high-pt tails of 
hadronic single particle spectra ("jet quenching") [ ||, ^] 
is thus a tentative signal for the formation of a decon- 
fined state. This possibility has motivated many stud- 
ies of medium-induced gluon radiation in recent years [ 

i, i, @, i, i, 0, m m ei i@@ @. 

Most recent studies of the non-abelian energy loss are 



carried out within the Gyulassy-Wang (GW) model [ ||] 
which mimics the medium by a set of coloured static 
scattering centers. Despite its simplicity, this model is 
of phenomenological interest, since the medium- induced 
radiative energy loss belongs to those observables which 
depend mainly on the average transverse colour field 
strength encountered by the parton rather than on 
the model-specific details with which this colour field 
strength is described [ [l9| . 

First studies of the GW model [ [l], |^| focussed on 
the rescattering of the hard quark and arrived at a 
radiative energy loss dE/dx — const independent of 
the path length. Baier-Dokshitzer-Mueller-Peigne-Schiff 
(BDMPS) [ |(| established later that gluon rescattering 
diagrams give the dominant contribution. They found 
a potentially dramatic linear increase dE/dx oc L of 
the energy loss with the medium thickness L. This 
can be understood in terms of an uncertainty argument 
of Brodsky and Hoyer [ ^| which relates the average 
transverse gluon momentum (k^) to the radiative en- 
ergy loss dE/dx oc (k^). Brownian motion of the rescat- 
tering gluon then implies (k^_) oc L and the quadratic 
L-dependence of radiative energy loss in the BDMPS for- 
malism. While all these studies start from the complete 
set of multiple scattering diagrams, Zakharov [@, || has 
advocated a different and very elegant approach to the 
same problem. In his path-integral formalism, the radia- 
tion cross section is determined by a dipole cross section 
which essentially measures the difference between elastic 
scattering amplitudes of different projectile Fock state 
components as a function of impact parameter. Baier, 
Dokshitzer, Mueller and Schiff (BDMS) have shown [0 
that the evolution of the rescattering amplitude in the 
BDMPS-formalism is determined by Zakharov's dipole 
cross section. 

The above mentioned calculations do not have a clear 
connection to experiment since the total radiative energy 
loss dE/dx is not an experimental observable. Unlike 
the situation in QED where the charged projectile can in 
principle be arbitrarily well separated from its radiation, 
medium-induced QCD bremsstrahlung is only an observ- 
able to the extent to which it is emitted kinematically 
well separated outside the typical hadronization cone of 
the hard parton. Realistic energy loss estimates thus 
require knowledge about the kx-differential gluon radia- 
tion spectrum. First calculations of this observable were 
published a year ago in the Zakharov [ and in the 

BDMPS-formalism [ [l4|. These studies contain essential 
steps towards a phenomenological application: especially, 
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Ref. [ |l2j shows how to rewrite Zakharov's path-integral 
formalism in a simple, numerically accessible form, and 
Ref. [ [l4| gives results for the radiation outside the kine- 
matically unresolvable hadronization cone of the parton. 

However, calculating the angular distribution is not 
the only problem in making contact with phenomenol- 
ogy. Another problem is that the concept of a homo- 
geneous medium of finite extent, underlying the BDMPS 
and Zakharov formalisms, may not be applicable to heavy 
ion collisions at RHIC and LHC. Unlike the situation 
in QED, where even the thinnest targets probed in ex- 
periments of the Landau-Pomeranchuk- Migdal effect 
amount to » 10000 — 100000 small-angle scatterings, one 
expects for the medium-induced radiation off hard par- 
tons in relativistic heavy ion collisions a much smaller av- 
erage number of rescatterings, say « 2 — 10 [ [lTj . A com- 
parison of the BDMPS- and Zakharov- formalisms with 
results for a small fixed number N = 1, 2, 3, ... of rescat- 
terings is needed to understand to what extent the con- 
cept of a homogeneous medium can be justified for such 
an extremely thin medium. The only existing studies [ 
p"2| , p"7| ] of scenarios with very few A < 3 scattering cen- 
ters assume implicitly an exclusive measurement of jet, 
gluon and recoil target partons. They arc thus restricted 
to the study of a subclass of all available rescattering di- 
agrams, and their results do not compare directly to the 
BDMPS and Zakharov formalisms (see below for more 
details). It is one of the results of the present paper to 
calculate the gluon radiation spectrum for a "medium" of 
very few N < 3 scattering centers, including all available 
rescattering diagrams, thus allowing for a direct compar- 
ison. More generally, we discuss in which sense the A-th 
order term in the opacity expansion of Zakharov's result 
is related to radiation off a target of a fixed number A 
of scattering centers. 

From a technical point of view, the present work is 
an application and extension of methods developped in [ 
|l9| |. There, we have derived the non-abelian Furry ap- 
proximation for the wavefunction of a hard parton which 
undergoes multiple rescattering in a spatially extended 
colour field. This Furry wavefunction provides a com- 
pact shorthand for the high-energy limit of the complete 
class of final state rescattering diagrams. Relevant for 
the present work is that we have developped in Ref. [ 
|T9| diagrammatic tools to determine under which condi- 
tions observables of a multiple partonic rescattering pro- 
cess are colour trivial. Colour triviality is the remark- 
able fact that for some rescattering processes, the sum 
of all contributions to the AT- fold rescattering depends 
on a unique A-th power of the SU (3) Casimir opera- 
tors, rather than to depend on colour interference terms 
associated to more than one colour trace. This renders 
the calculational problem essentially abelian and is cru- 
cial for going beyond very few rescatterings where brute 
force perturbative calculations are still feasible. Colour 
triviality is thus an important property in the study of 



non-abelian gluon bremsstrahlung. 

The present work applies the non-abelian Furry ap- 
proximation to the calculation of the medium-induced 
gluon radiation cross section. In section ||, we derive 
the radiation cross section in terms of Furry wavefunc- 
tions. From this, we derive in section |l| a set of dia- 
grammatic identities which automate the proof of colour 
triviality. In section IV, these identities are applied to 
derive an explicitly colour trivial diagrammatic expres- 
sion for the A-th order term of the opacity expansion 
of the k_L-differential gluon radiation spectrum. Contact 
with analytical expressions is made by showing that the 
sum over all A-th order contributions results in a path- 
integral expression which essentially coincides with Za- 
kharov's result. In section [v|, we relate this path-integral 
expression to the radiation off a target of a fixed num- 
ber of scattering centers. All calculations in sections |l]- 
^ are for gluon radiation off an idealized free incoming 
quark which satisfies plane wave boundary conditions at 
far backward position. The corresponding radiation cross 
sections are denoted by a superscript (in) . In section VI , 



we extend these calculations to the radiation off nascent 
quarks produced in the medium. The gluon radiation 
cross section associated with them is characterized by 
a superscript (nas) and gives access to the interference 
pattern between hard and medium-dependent radiation. 
The main results and further perspectives are finally dis- 
cussed in the Conclusions. 



II. THE FORMALISM 

We consider gluon radiation off a hard quark which 
undergoes multiple elastic rescattering in a spatially ex- 
tended colour field. The gluon of momentum k = (u>,k) 
carries away a fraction x of the totally available energy 
of the initial parton, u) = xE\. A typical contribution 
to the radiation amplitude is shown in Fig. |l|. For the 
coloured potential describing the medium, we choose the 
ansatz of the Gyulassy-Wang model [ 0| 
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(2.1) 

(2.2) 
(2.3) 



The static scattering potentials ipf(x) show a sufficiently 
rapid (e.g. Yukawa-type) spatial fall-off. At the i-th 
scattering center, a specific colour charge d = di is ex- 
changed. T d , (d = 1, . . . , N% — 1), denotes the generators 
of the SU(N C ) fundamental representation, and the to- 
tally antisymmetric structure constants f a bc denote the 
adjoint representation. The latter appears in the rescat- 
tering amplitude of the emitted gluon. 
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FIG. 1. Contribution to the gluon radiation amplitude 
(2.4), involving M-, L-, and JV-fold rescattering of the in- 



and outgoing partons respectively. 

Summed over arbitrary many rescatterings, the gluon 
radiation amplitude M- a _ tbc in this extended colour po- 
tential takes the form 

M a ^ bc = J d 3 yI aai (y, Pl ) (-» S . 7|tl (T Cl ) aibl ) 

xe-^I^ c {y,k)h lb {y, P2 ). (2.4) 

Here, the term e -6 ' 2 " is the adiabatic switching off of 
the interaction term at large distances. In our calcula- 
tion, the limit e — > does not commute with the longi- 
tudinal ?/;-integral, and this regularization has to be car- 
ried through intermediate s teps of the calculation [ [l2| . 
To explain the notation of (2.4), we consider the typical 
diagrammatic contribution given in Fig. | The TV-fold 
rescattering of the quark, emitted from the radiation ver- 
tex with momentum P21 and colour b\, and appearing in 
the final state with momentum p2 and colour b is deter- 
mined by the component of [see Fig. [l] for details of 



notation] 



y,?2) 



x i(^ + m) 70 [ _. A ( g)(xj)] 



lXl .(p 2 , I+1 - P2 , I ) N ) w(p a ). (2.5) 



and the rescattering effects on the emitted gluon are 
taken into account by the component Ic-lc of 



xH4 5) (x t )l 



-i g»* ^ 



kf — m 2 + i e 



x e -ix<-(ki + i-k<n e ^ L 



(2.7) 



where e^ L is the polarization of the gluon. To calculate 
the high-energy limit of the radiation spectrum, one has 
to keep on the amplitude level the leading order in norm 
and the next-to-leading order in the phase. In this limit, 
other scattering contributions (e.g. those involving the 
4-gluon- vertex) are negligible [ and the expressions 
(2.5)- take the form of non-abelian Furry wavefunc- 
tions [ |19[. In appendix we show that by approxi- 
mating 7 M1 ^(y, k) to leading order 0(\/u) in the norm 
and next to leading order in the phase, the correspond- 
ing sum over arbitrary many rescattering centers takes a 
particularly simple form 



i^( y ,k) = J2i ftl(L) (y,k) 

L=0 

= J dx ± G (9) (y;x|^)F(x,k). (2. 



Here, F(x, k) denotes the asymptotic plane wave 



F(xj_,a:;,p2) = cxp \ -ip^ • xj_ + i x< 



2p 2 



(2.9) 



and the Green's function G( g \ is an approximate solu- 
tion of the non-abelian Dirac equation in the spatially 
extended colour potential A^K This Green's function 
describes the leading transverse deviation of the rescat- 
tering parton from a straight line propagation. It can 
be represented in terms of a path-integral over a path- 
ordered Wilson line Wi g \ ([r]; Z, z'\ , 



Here, the path-ordering V implies that Aq (x i+ i ) stands G( g ) (r, z; r', z'\p) = 



to the right of A^f 1 (jq ) , and the momentum transfers to 
the quark line are written as Fourier transforms of the 
static scattering potential with respect to the relative 
momenta P2.i+i— P2,i- «(P2) i s the spinor of the outgoing 
quark. In complete analogy, rescattering effects of the 



incoming quark are described by the component I a 

i (M) (y,Pi) = vH Pl )v n 



(M) 



of 



pi 



(2tt) 3 

>i,i + m) 70 
— m 2 + u 



VriO rx,, <( | / d£r 2 (£) } W([r];z,z') , (2.10) 



W (g) ([r};z,z')=Vexpl-i J d£A^(r(0,0 1 • (2.11) 



In what follows, we shall explicitly denote or supress the 
longitudinal coordinates of these Green's functions, de- 
pending on whether this information can be easily in- 
ferred. 



x e 



•(Pl,i+1— Pi,-, 



'Pi,i-y 



7o : 



(2.6) 



In analogy to ( 2J3 ) , we have shown in a previous paper [ 
|l9| that the sum over an arbitrary number of final state 
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rescatterings of the outgoing quark takes the compact 
form 



l(y,p 2 ) = £ i {N) (y,P2) 



N=0 

= e- l ^y'D 2 J dx ± G {q) (y; X \p 2 )F(x,p 2 )v(p 2 ) . (2.12) 

As denoted by the subscript (q), the Green's function 
contains in this case a Wilson line in the fundamental 
quark representatio n of SU(N C ). The only new ingre- 
dient compared to (2.8) is the differential operator L>i [ 
0,0 



Di = \-i 



a-V a ■ (pi - npi) 



2 Ei 2 Ei 

cc = 7o 7 ; z = n ■ x 



Pi 



(2.13) 



This operator contains in a configuration space formula- 
tion the 0(1/ E) corrections to the propagator of momen- 
tum p2i- Inclusion of these corrections for the fermion 
propagators entering the radiation vertex is necessary, 
since the leading l/i?-contribution cancels due to the 
Bloch-Nordsieck structure of the radiation vertex [ [l!| . 
In terms of the Green's functions G/ q -\ and G( 9 ), the 



gluon radiation amplitude (2.4) takes the form 



The radiation probability (2.17) has a simple diagram- 
matic representation in configuration space which we give 
in Fig. ^. We denote by solid lines the full quark Greens's 
functions G( 9 ) , and by dashed lines the gluon Green's 
functions Gi g y Contributions to the amplitude are on 
the l.h.s., those to the complex conjugate amplitude on 
the r.h.s. . We consider only the case that the gluon ra- 
diation occurs at larger longitudinal position in the com- 
plex conjugate amplitude, y > y . Taking twice the real 
part in ( 2.17 ) accounts automatically for the other case. 
Transverse integration variables and colours are speci- 
fied at the end points of the Green's functions. In this 
way, the simple diagram of Fig. || characterizes the rather 
lengthy expression (2.17) completely. 
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FIG. 2. Dia grammatic representation of the radiation 
probability (2.17) in configuration space. For details, see text. 



M a . 



dyi e % 



\vi\ 



d 2 y 



xG^(x i; y| Pl )r y (T c *) ol6l Gg) c (y;x a |u;) 



( ^(y;x 2 |p 2 )e 2 " 2 « ^ 



(2.14) 



Here, we have used the fraction x of the incident energy 
carried by the emitted gluon, w = x E± , to write 



The following sections are devoted to a study of the 
k_L-differential gluon radiation cross section off a free in- 
coming quark. From ( [2.17 ), we obtain 



d 3 a (in) 



d{\nx)dk_ L (2tt) 2 4£?(1 - x) 2 



Cp rc 2 Re 



J dpt(\M a ^ bc \ 2 ) 



dyi / dyi e 
' Vl 



iq{yi-yi)-e \vi\-* \vi\ 



dy dy dxi dx.\ dx„ dx g dx 2 e 



ik±-(x g — x 9 ) 



x m„ 



q = Pi - P2 - W 



2(1 -x) E l 



(2.15) 



x(G-(x 1 ;yb 1 )T- fci [-G^(y; Xg \ 



Also, we have introduced the interaction vertex f as the 
notational shorthand for the spinor structure of the am- 



plitude ( |2.14| ), 

f y =v^[px)Dx'y e-"fD2v(p 2 ). 
This leads to the radiation probability 

(\M a ^ bc \ 2 ) = 2Re [ + d Vl [ + dy ie ^y-^ 



-« \yi\~ e \vi\ 



(2.16) 



dy dy dxi dx± dx 2 dx 2 dx g dx g T y f y 



xe 
xC 

xG-5(x 2 ; y | P2 )G-(x s ;y| W )7- ai G^(y:x, „ 



xG (,°) (xr,y|pi)T^ 6i G^ } c (y; x g \u) G b {q) (y; x 2 |p 2 ) 



xe 



ipi_L-(xi— Si)- ip2j_-(x 2 — x 2 )— ikj_-(x g — x 3 ) 



xG^(y;x 2 |p 2 )G^(x 2 ;yb 2 ) 
d 



x ( -G-(x 3 ;y|^) ) G^{fMPi) ) , (2-18) 



where 



Cprc — 



(2tt) 2 x 2 E 2 (1 - x) 2 



(2.19) 



accounts for a com bina tion of recurring kinematical fac- 
tors. Expression ( j2.18j ) is given in the frame in which 
the incoming quark has vanishing transverse momentum 
Pi_i_ = 0. The medium average (. . .) over the colours 
and posit ions of the scattering centers is specified in sec- 
tion II A below. In appendix pi we give details of why 



(2.17) the combination of interaction vertices (2.16) can be writ 
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ten to leading order in 0{x) as derivatives acting on the 
gluon Green's functions. 



A. Medium average 

The main tool of our analysis of the radiation cross sec- 
tion 02.180 wil1 be 

an expansion of the Green's functions 
in powers of the scattering potential Aq: 



z 

G(r, z; r', z'\p) = G (r, z; r', z'\p) - i J d£ 

Z 

x J dpG {r,z;p,£,\p) A (p,£) G {p,&r',z'\p) 

XL XL 

+V l ^ l ^ l dPldp 2 G °( r ' 2;;Pl '^ 1 ^ 

x i A (px , £1) G (Pi, fi ; Pi, 6 \v) 
xiA (p 2 ,&)G(p 2 ,&;r',z'\p). (2.20) 

Here, Gq is the free non-interacting Green's function 



G (r,z;r',z'\p) 



P 



2iri(z' — z) 



exp ■ 



ip (r 



2{z' - z) 



(2.21) 



and the only distinction between the quark and gluon 
Green's function lies in the colour representation of the 
spatially extended colour po tentia l ^4o- An expansion of 
the radiation cross section ( 2.18 ) to order 0(Aq u ) will 
turn out to be an expansion to n-th order in the opacity 
of the medium [ [l2], |9). Accordingly, we shall refer to 
( 2.20 ) as opacity expansion. 

We no w discuss the calculation of the medium average 
in fl2.1§| ) to fixed order in opacity. To this end, we write 
for a single scattering potential, centered at (fj,2j), 



xe 



(27p 



(2.22) 



Here, the high-energy approximation &o(<7) « ao(q±) im- 
plies that the effective momentum transfer occurs at fixed 
longitudinal position £ = z- L . The medium average (...) 
is the average over the transverse and longitudinal po- 
sitions (ri,Zi) of the scattering potentials and over the 
colours cij exchanged with the i-th scattering center, 



(/> = 



1 



N 



\i=l cii 

x/(fi, ...,r N ;zx,.. . . .,a N ) . (2.23) 



A±_ is a total transverse area which we divide out to 
regain the cross section per unit transverse area. N is 



the number of different single scattering potentials up to 
which the function / is expanded. 

To simplify notation, we re place in what follows the 
discrete sum over Zi in ( 2.1C ) by an integral over the 
density n of scattering centers 

Ao(x x ,£) = J dzin(zi)tp?(x ± ,OT a . (2.24) 

Since the effective momentum transfer from the single 
potential iff (xj_ , £) occurs at £ = zi , wc shall often work 
with £ as integration variable. 



B. Leading O(x) approximation 



We restrict our analysis of ( 2.18Q to the kinematical 
region where the fraction x of the energy carried away 
by the gluon is small, x <C 1. This 0{x) approximation 
is commonly used in the analysis of the abclian and non- 
abelian LPM-effect. It focusses on the kinematical region 
in which most of the radiation occurs. Technically, it 
allows for two important simplifications: 

1. The leading 0{x) spin- and helicity-averaged prod- 
uct of the interaction vertices takes in momentum 
space the simple form 



1 y 1 y 



— ku_ ■ ku_ 



(2.25) 



where ku_ and ki^ are the transverse momenta of 
the gluon at the radiation vertex in the amplitude 
(ku_) and complex conjugate amplitude (ku_). In 
configuration space, these momenta are conjugate 
to the y-derivatives of the gluon Green's functions. 
Details are explained in appendix |b|. 

2. Propagation of the transverse plane waves by 



Green's functions leads in ( 2.17 ) to the appearance 
of longitudinal phase factors: 



J dx 2 G (xi,zi;x 2 ,Z2\E)e , pj -' 



-ip_L-Xi 



exp 



-t — — [Zo — Zl 



(2.26) 



To leading O(x), the longitudinal phase receives 
only contributions from the transverse energy 
p^/2w of the gluon. Transverse energies p 2 ± /2Ei 
of the quarks are suppressed by one order in x. In 
this sense, the leading 0(x) approximation neglects 
the Brownian motion associated with the rescatter- 
ing of the hard quark. Gluon rescattering is the 
dominant contribution, as observed by BDMPS. 
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III. OPACITY EXPANSION 

For more than N = 1 scattering centers, the opacity 
expansion of the radiation cross section (2.18) involves 
non-vanishing interference terms between amplitudes of 
different powers in Aq. Two examples of non- vanishing 
interference effects between amplitudes of order O(Aq) 
and O(Aq) in the N = 2 opacity expansion are e.g. 



® © © 
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def 



def 



- x x- - 



(3.1) 



In what follows, we use the expression contact term to 
denote scattering centers which link with two gluons 
to one amplitude [as e.g. the first scattering center in 
( |3 .![ )]■ The notion real term characterizes a scattering 
center which exchanges one gluon with the amplitude 
and one with the complex conjugate amplitude [as does 
the second scatterin g cen ter in (3.1)]. The structure of 
the medium average ( |2.23 ) ensures that all diagrammatic 
contributions are combinations of real and contact terms. 
For a medium of density n(£) and thickness L, arbitrary 
combinations of N real and contact terms contribute to 
the same order (a s n(£) <i£) in opacity. Hence, con- 
tact term contributions cannot be neglected in the calcu- 
lation of ( 2.18p . Calculations which do not include con- 
tact terms [jl2|, [lTj describe an exclusive measurement in 
which the number of recoil target partons is determined. 
In section [II A, we discuss the fundamental properties 
of contact terms as well as a set of identities which re- 
late them to real terms. In section 1IIB. we use these 



identities to derive a manifestly colour trivial expression 
(3.13) for the radiation cross section off N scattering cen- 
ters. Then, we calculate (3.13) in an opacity expansion 
for N = 1 and N = 2. 



A. Notation and identities 



The medium average ( 2.23| ) ensures that contact terms 
do not transfer a net momentum to the final state. 
However, they can change the transverse momentum 
between the outgoing particles [ |l^]. This can be 
seen, e.g., by writing for the first scattering center in 



(3.1) the corresponding cont ribut ions ( 2.22 ) and realiz- 
ing that the medium average (2.22) results in a ^-function 
<^ 2 ' ) (<Zi_l + with the transverse momenta of both 
gluons appearing in the same amplitude. 

Contact terms satisfy several important identities. All 
the following identities can be established by writing out 



the opacity expansion ( 2.20| ) for the last scattering cen- 
ter before the cut and employing the medium average 
( p. 23 ). Several examples of this technique were discussed 
explicitly in [ [HJ. Here, we consider diagrams to N-th 
order in opacity. We denote by a shaded box the longi- 
tudinal region in which the N — 1 first interactions take 
place, and we specify the contribution of the last N-th 
scattering center explicitly. We start with the case that 
the last interaction links to the quark line. Irrespective 
of whether it occurs after yi , 



(3.2) 



or before j/;, 



(3.3) 

the real contribution equals in both cases minus the two 

corresponding contact terms and thus cancels against 

them. We note that (3.2) holds exactly. It is based on the 

identity derived in Fig. 7a of Ref. [ Equation (fr.Sj), 

in contrast, involves the leading 0(x) approximation: the 

r- — . To 2 1 

l.h.s. of (3.3) involves a phase exp —ij^{yi — yi) whose 

transverse momentum differs by the amount of 
the last rescattering from the corresponding phase on the 
r.h.s. of (|3.3[). This phase does not contribute in the lead- 



ing 0(x) approximation where the momentum transfer to 
the quark line is neglected and only transverse energies 
of the gluon are kept. The same argument allows for an 
identity concerning the first scattering center: 



(3.4) 

The important implication of this identity is that the sum 
of all rescattering contributions to (2.18) which have the 
first interaction before yi, vanishes. This eliminates a 
large class of diagrams from our calculation. 

In addition, there are identities in which one of the 
gluon lines is touched: 



G 



-X- - 



JC- - 



(3.5) 

These identities are exact and were derived already in 
Fig. 7b of Ref. [g. 

So far, we have discussed identities for which the last 
interaction occured either after yi or before yi. Now we 
give relations which hold for a last interaction at longi- 
tudinal position yi < £ < y\. These identities necessarily 
involve knowledge of the colour algebra. The simplest 
example is 



(3.6) 

Here, the r.h.s. shows a star over the interaction term. 
This notation indicates an interaction term without 
colour factor. The colour information is absorbed in the 



prefactor of the diagram. To derive (3.6), we denote the 
colour generators of the left and right shaded blocks in 
( |3.6[ ) by factors L and i?, respectively. The colour trace 
of the first term on the l.h.s. then takes the form 



Tr [LT d T c T d R] 



C, 



Ca 
2 



Tr [L T c R] , (3.7) 



where T d is the generator associated with the N-th inter- 
action and T c that of the quark-gluon emission vertex at 
yi. The contact terms on the left and right hand side of 



(3.6) are proportional to Tr [L T c R] with prefactors 



2 



and — %S respectively. In the leading 0(x) approxima- 
tion, there is no complication from the momentum struc- 



ture of the diagrams and equation (3.7) thus ensures the 
identity (3.6). 

In a similar way, one can establish an identity for dia- 
grams in which one of the interactions links to the gluon 
line: 



I X 



- x 



2^A 



This is a direct consequence of the relation 
i fede Tr [L T c T d R] — ^ Ca Tr [L T s R] 



(3.8) 



(3.9) 



B. Colour triviality of the gluon radiation cross 
section 

The identities of the last subsection have two impor- 
tant properties which we exploit systematically in what 
follows: (i) they allow for many diagrammatic cancella- 
tions between real and contact terms, and (ii) they ensure 
that the remaining terms combine to expressions which 
are proportional to Cj[- 



1. Classification of diagrams 

In what follows, we consider all diagrammatic contri- 
butions which are N-th. order in opacity. The N — 1 
first interaction terms are again summarized diagram- 
matically in a shaded block on both sides of the cut. 
The last interaction is specified explicitly. We consider 
three cases 

Case I : The last interaction occurs at longitudinal po- 
sition £at < yi . 

In this case, the identity ( |3.2| ) ensures that the sum of all 
real and contact terms of the last interaction vanish. 

Case II : The last interaction occurs at longitudinal po- 
sition y t < £jy < Vi- 
la, this case, we have two real and four contact contribu- 
tions from the last interaction. According to the identi- 
ties (3.6) and (|3.8|), they add up in the following way: 



- -x 



(3.10) 

In the last line, we have introduced a simple short- 
hand which summarizes all non-vanishing contributions 
for yi < £jv < Vi- F° r what follows, it is important that 
this shorthand comes with a unique colour prefactor Ca 
which absorbs the colour algebra of the last interaction. 

Case III : The last interaction occurs at longitudinal 
position yi < £ N . 

In this case, we have four real and six contact term con- 
tributions. Identity ( |3.2| ) ensures that the one real and 
two contact terms shown there cancel each other. Iden- 
tity J3.5|) leads to the cancellation of two real versus two 
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contact terms. For the remaining one real and two con- 
tact terms, we introduce the notational shorthand: 



(3.11) 

Again, the sum of all non-vanishing contributions comes 
with a colour prefactor C 'a which absorbs the colour al- 
gebra of the last interaction. 



2. Iteration 

The above classification of diagrams shows that the 
colour structure of all surviving contributions to the last 
iV-th interaction can be absorbed in a prefactor Ca- This 
allows for an iteration of the above procedure to the ear- 
lier (N - l)-th, (N - 2)-th, (N - 3)-th etc. interactions. 

We consider first the case that the m last interactions 
occur at positions larger y h i.e., yi < £ N - m < CiV-m+i < 
... < ^at. For the J V-th interaction, we are left with 
the diagrams of ( 3.1l| ). Since the colour structure is ab- 
sorbed in the prefactor Ca, the arguments leading to 
( |3.11 ) apply also to the (N — l)-th interaction. If the 
(N — (m + l))-th interaction also occurs after yi, then 
the sum of all nonvanishing contributions is given by the 
recursively defined diagram 





m+1 











(3.12) 

For iV-fold rescattering, an arbitrary number m < N in- 
teractions will occur after y\ , the remaining (N — m) in- 
teractions will occur at longitudinal positions between y\ 
and yi. The contributions of all possible diagrams with 
interactions before y\ cancels due to the identity (3.4). 
The radiation cross section ( 2.1S| ) to N-th. order in opac- 
ity can thus be represented in the compact diagrammatic 
and explicitly colour trivial form 



dV m )(iV) _ a s 
d(\nx) dk ± ~ (2tt) ; 



—z Nc Cf C a 



N 



(3.13) 



N 

I" 



This is the diagrammatic manifestation of the observa- 
tion of BDMPS that gluon rescattering alone gives the 
l eadin g 0(x) contribution t o the radia tion cross section 
(|2.18|) . As can be seen from (ft.lCD and (|3.12|) , all N scat- 
tering centers involved transfer momentum to the gluon 
line. 



We introduce for the diagrammatic part of (3.13) an 
analytical shorthand 



I' 



Y «'[gi,£»]ig[i,iv] 



( ^G(y;x 9 |c; 



dy dy ds dx g dx g 

G (y;y> 2 ) 



xe 



ik±-x a 



— G(x s ;y|w)] CVi y.y,: s. i . (3.U) 



To specify our notation, we com pare ( 3.13 ) and ( 3.14 ) to 

l 



the full radiation cross section (|2.18|) . The term 



in ( [3.12 ) is the leading order 0(x) of the prefactor C wc . 
The colour algebra of (2.18) simplifies to the Casimir fac- 
tors in (3.13). Accordingly, the functions G introduced 
in (3.14) are abelian. Some quark Green's functions in 
the diagram of ( [3.14 ) are non- interacting to arbitrary or- 
ders in opacity. This allows to integrate over Xi , xx and 
x 2 in ( p. 18 ). The new integration variable s in (3.14) 
denotes the transverse coordinate of the quark at longi- 
tudinal position yi in the complex conjugate amplitude. 



The longitudinal yi- and yj-integrals in (2.18) are given 
by the notational shorthand 



/ = 2Re 



dyi 



dfre-'M-'Mf, (3.15) 



where the function / denotes an unspecified integrand. 
To N -th order in opacity, the medium average (...) in 
fl2.18|) results in N transverse momentum q^-integrals, 



f 



(2ny 



\ao{qi±)\ f, 



(3.16) 



which average over the elastic scattering cross sections 
\ a o(Qi±)\ 2 f° r eacn scattering center. Also, it leads to 
N integrals over the allowed longitudinal positions £j G 
Kr m 'Cr ax ] of these N scattering centers, 



f = 



<%i n (t,i) f ■ 



(3.17) 



These integrations are combined in (3.14) into 



JV 

/ /= n/ / / 



(3.18) 



s 



As a consequence of factorizing in ( p. 18 ) the colour al- 
gebra, the elastic cross sections \ao(qi±)\ 2 and the lon- 
gitudinal momentum integrals ov er the only rem- 
nants of the opacit y exp ansion (|2.20| ) which are left 
in the integrand of ( |3.14| ) are transverse phase fa ctors . 
Hence, the functions G are defined by replacing in ( 2.20| ) 
Aq(p, £i) — > e~ lqi ± p . For each qi±, two phases appear in 
the integrand of ( 3.14 ), one being the co mplex conjugate 
of the other. To gain familiarity with ( pM4| ), one may 
consider simple examples like 



J J J ^il(Ci) J dy dy ds dx g dx g 
d 



9y 



J dpi Go (y; pi , £1 | w) e"^ ^ G ( Pl , & ; x fl |«) 



Ce -*x-(x s -x 9 ) Go(y .-| p2) I ii Go (x ff ;y|, 







x / dpiGo(y;pi,abi)e l9li ^ Go (pi, &; s, j«|pi) , 
= / / (k ± +qix)-k± / dan(6) 

JYJqx J yi 



Xe -i^(vi-^)-i {kX+ 2l 1X) (falft) 



(3.19) 



Here, we have used again the leading O(x) approximation 
to neglect the subleading contributions to the longitudi- 
nal phase of the integrand. 

All N-ih order contributions to the radiation spectrum 



(3.13) can be written as a sum of products of the form 
dV m )(A) a s 1 



A' 



d(lnx)dk x (2tt) 2 uj 2 



N C C F C% 



J qi J qn 



Here, the factor 



Z(i) - 



(3.20) 



(3.21) 



denotes an integral over longitudinal phase factors. To 
leading 0(x), this phase is obtained by iterative use of 



(2.26). The factor in (3.21) is the integral of this 



phase over the allowed longitudinal position of scattering 
centers times a factor —1/2 for eac h con tact term linking 
to the gluon line. In the example ( 3.19 ), is given by 
the £i-integral. 



The factors i n (3.20 ) denote the results of the par- 
tial derivatives in ( 3.14 ) for the i-th diagrammatic con- 
tribution. They are easily read off from the diagrams by 
adding up the momentum transfers to the gluon line in 
both amplitudes. For the example ( 3.19 ), this leads to 
the fact or (k j_ + Qix) ■ kj_. We now turn to an evalua- 



N = 2. In section IV, we derive then an all-order expres- 
sion ( 1.15| ) for this radiation cross section. As explained 
in sections [v] and VI, this all order expression gives a 
much simpler access to the opacity expansion since it in- 
volves only N + 1 terms to iV-th order. 



C. N = 1 Gunion-Bertsch radiation spectrum 

Here, we derive the Gunion-Bertsch radiation spec- 
trum [ |22| as the N = 1 term in the opacity expansion 



of ( 3.13| ). This provides a first consistency check for our 
formalism and allows for a simple illustration of our di- 
agrammatic rules. The N — 1-term of (3.13) has four 
contributions 



dV in )(iV = 1) a s 1 



d{lnx)dk ± (2tt) 2 uj 2 

x ( — : _ + 2 

+ — r* "" : t — - + 2 



Nc Cf Ca 



a 1 

pj^ja ~ji Nc Cf Ca H r « Z (i) 



(3.22) 



Here, we have labelled the different diagrammatic contri- 
butions by (1) - (4). Our diagrammatic rules distinguish 
between interactions before and after yi. For that rea- 
son, the diagrams (3) and (4) are both includ ed in ( 3.22 ). 
Contribution (2) is the example studied in (3.19). The 
T-factors are read off from the i-th diagramm by adding 
up the transverse momentum flows into the gluon line on 
both sides of the cut: 



r ( i) = (kx + q 1± ) ■ (k x + <?i_l) 
r (2) = (k_L + q 1± ) ■ k± 



r 



(2) 
(3) 



(4) 



(3.23) 
(3.24) 
(3.25) 



For th e leading 0(x) contributions to the factors Zu\ in 
( |3. 19 ) , only the transverse momenta of the gluon have to 
be considered. In ( 3.22) ) , the gluon line comes with two 
possible transverse energies 



k 2 



Qi 



(k ± + q u 

2lo 



(3.26) 



Including a factor —1/2 for each contact term linking 
twice to the gluon line, we find for the phase factors 



tion of (3.12) in the low opacity expansion for N = 1 and 



$(3)(yi»j«) 
®(4)(yi,Vi) 



' dtn(£)e- zQ{y '- yi) , 
dtn(Z)e- iQ ( yt - y,) . 



(3.27) 
(3.28) 
(3.29) 
(3.30) 
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For a homogeneous density distribution = no of a 
medium of thickness L, the corresponding factors Zn\ 
given in ( 3.21) are 



Z (i) - 

Z (3) = 



Qi ' 

2 L no 



(4) 



-2 L no 

QQi 

-Ln 



(3.31) 



Inserting this into the last line of ( 3 . 22| ) , we find 

dV™)(iv = l) 



G?(ln x) dk j 



(Lno) N C C F C A 



n 2 



k 2 ± (qi± + k±_Y 



(3.32) 



This is the leading 0(x) contribution to the Gunion- 
Bertsch [ ^2| radiation cross section times the opacity 
Lno °f a homogeneous medium of thickness L. 



D. N = 2 Opacity expansion 

Before turning in the next section to a systematic all- 
order analysis of the radiation cross section (3.13), we 
calculate here the N — 2 contribution to the opacity 
expansion. 



+ 2- 
+ 

+ 2- 
+ 

+ 2- 





-x- x- 


X Xr 




(1) 




-b- x- 


* 

x - - 




(2b) 


,- x - 


- - x- 


x - - 




(4) 








— x 




,-e- - 


- - x- 


X - - 




(2c) 


,x-x 








(6) 








X x — 




,-e-x- 






X 


(5c) 



+ 2- 
+ 4 
+ 2- 
+ 4 
+ 2- 
+ 4 



X Q- 



(3b) 



(5b) 



FIG. 3. All non-vanishing con tribu tions to the N — 2 radi- 
ation cross section as given by (3.13). Diagrams are labelled 



with a double index (ix), see text for details. 



All terms contributing to the N — 2 contribution of 
(3.13) are listed in Fig. |j. We have labelled them with a 
double index (ix). The first index i = 1, 2, . . . , 6 indicates 
how the diagram contributes to the interaction vertex 

2 



r ( i) = 

r (2) = 
r (3) 

T(4) = 

r (5) = 

r (6) = 



(ki 
(kj 
ki 
(kj 
(kj 
(kj 



<?i_l + Q2j 
Qix) , 



9i_L + <72±) ■ (kj 

qi±) ■ k_L , 

<Zl_L + <?2x) • k_L 



Qi±) 



(3.33) 
(3.34) 
(3.35) 
(3.36) 
(3.37) 
(3.38) 



The integration variables q±± and q>2i_ for the first and 
second momentum transfer can be relabelled freely. As 
a consequence, we have to consider only three different 
transverse gluon energies 



k 2 



ii = 



(k ± + q 1± ) 2 



2to 



(kj, + gl_L + g2jj 

2w 



(3.39) 



and no t erm o c (kj _ + <72_l) appears in the interaction 
vertices (^33|)-(^38|). 

The phase factors can be read off from the diagrams of 
Fig. H as described above. They are now integrals over 
the allowed positions £i and £2 of both scattering centers, 
e-g-, 



xe 



'VI 

-iQ(v 



-6) 



d6n(a)»(6) 

i(fa-fi)-»Q2 (£1-1/1) _ (3.40) 



The second index x — a, b, . . . labels in Fig. [3] the differ- 
ent diagrammatic contributions to the same term Tuy 
Denoting by the phases the sum over all contribu- 



£36) 



tions (e.g. $ (3) = $ (3o) 
factors Z(j} according to ( 3.21 ) 

z L " 

Zr 



^(3c))i we nn d for the 



2Q 2 



-(2) 



1(3) 



-2 + 2 cos(iQi) 



Z (4) = 



Z (5) = 



-((i) 



-L 2 2 
2Q 2 "0, 

2-2 cos(LQi) - L 2 Q 2 
Q? Q2 

2-2 cos(LQi) + £ 2 Q 2 
-2 + 2 cos(LQi) 2 





(3.41) 




(3.42) 




(3.43) 


2 

> 


(3.44) 


2 

> 


(3.45) 




(3.46) 



These terms are multiplied by a factor —1/2 for each con- 
tact term linking twice to th e glu on line. The correspond- 
ing radiation cross section ( 3.20 ) shows L-dependent os- 
cillations which are due to the interference between the 
gluon radiation of the two spatially separated scattering 
centers. We consider especially for fixed opacity L no the 
limits L — > and L — > 00 in which the two scattering 
centers sit on top of each other and at arbitrary large 
separation, respectively. We find 



lim 

L — >oo 



e?(ln x) dk^ 

[R(k± - 



7T* 



N C C F C\ 



(Ln y 



L no —fixed 

qiX,q2±)-R{k±, qu _)} , (3.47) 
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and 



d 3 a( in HN = 2) 

lim 



" s at r r>2 ( Ln ») 
-zJ\ c C F C A — - — 



L rig— fixed 



L^o d(ln x) dkj_ 

x f f [R(k±,qi±+q 2 ±) -2R(k±,qi±)] . (3.48) 

We have written these limiting cases in terms of the 
momentum-dependent part R of the Gunion-Bertsch 



QCD radiation spectrum (3.22), 
RQt±,q±) 



kl (kj 



(3.49) 



which is associated to a quark of transverse momentum 
k^, receiving a momentum transfer q± in the scattering. 

One may expect that in the coherent factorization limit 
L — > 0, the radiation spectrum corresponds to one sin- 
gle Gunion-Bertsch term R(k±, qi± + q 2 ±) with effective 
momentum transfer (qix + q 2 ±), and that in the inco- 
herent limit L — > oo (which for QCD is of course only 
of formal interest), the radiation cross section is the sum 
of two independent Gunion-Bertsch terms. The results 



(3.48) and (3.49) differ from these expectations by the 
term 



-^jy c c F c A — - — 



2i?(k ± ,q lx ). (3.50) 



The origin of this term will be explained in section [v| 



IV. A PATH-INTEGRAL FOR ARBITRARY 
ORDERS IN OPACITY 

In this section, we derive a path-integral representation 
for the radiation cross section (3.13) which is valid to ar- 
bitrary orders in opacity. We start from (3.14), summed 
over powers of opacity: 



I' 



— - — / / dydydsdu ( -^-F(u — y' 
A± J Y J * y \dy K y ' 



d 



dr 3 (yi) 



M(r 1 ,T 2 ,r 3 \yi,yi) 



(4.1) 



We have separated in this expression scattering contri- 
butions from longitudinal positions before and after yi. 
The integration variables u and y denote the transverse 
positions of the gluon at yi in the amplitude and com- 
plex conjugate amplitude, respectively. All contributions 
in the region yi < £, < yi are given by 



oo „ 

M(r 1 ,r 2) r 3 |y / ,yO = E C '2 / 

n=0 •%<.£<]<g[l,' 



xG(r 3 (yi),yi;r 3 (yi),yi\u)) 
xGo{r 2 (yi),yi;r 2 (yi),y t \p 2 ) 
xG(ri(y ; ),yi;ri(yi)^/bi) 



(4.2) 



where the transverse paths rj(£), i E [1,3], satisfy the 
boundary conditions 



n(yO = s, 
r 2 (yi) = y , 

i"3(y/) = y , 



ri(y/) : 
r 2 (Vl) 



y, 
-y, 

- u . 



(4.3) 
(4.4) 
(4.5) 



The contributions coming from longitudinal positions £ > 
yi are summarized by the function 

F(U y) = f] C% f [d Xg d5i g 



xG(u,yr,x g ,z + \u) e 
xG(5c g ,z+;y,yi\u) . 



(4.6) 



Here, we have anticipated that the function F depends 
only on the difference u — y of the two integration vari- 
ables u and y. The consistency of our starting point (4.1) 
with the rescattering part to the radiation cross section 
(|3 .13) is straightforward to check by inserting ([h^)-([ll|) 
into ( fi~i"l) 



In what follows, we discuss how to simplify (4.2) and 
(4.6) in terms of the dipole cross section 



a(p) = 2C A J ^M^)! 2 (1-e-^ P) 



(4.7) 



1. Gluon rescattering for £ > yi 

We start with the analysis of the function F in ( |4.6| ) 
which describes the gluon rescattering at longitudinal po- 
sitions £ > yi. The first order term m — 1 reads 

C A 11 y 

= C A f_ + 'd£in(ei) J |^| ao ( 9± )| 2 

x J dx g d5t g e-^<^-^ /dr(Ci)dr(ei) 

x G (u, yi ; r (ft ) | w) G (r (ft ) ; x g , z + \ u) 
x G (x s , z + ; r (ft ) \u>) G (r (ft ) ; y , yi \ w) 
x ^"3( r (?i)-r(€i)) _ {J (4. 8 ) 

= -\ f_ + dft <ft) <r (u - y) e- lk -( u -^ , (4.9) 



11 



The factors in the last line of ( |4.8| ) correspond to the 
three terms on the right hand side of (|3.11 ) and can be 
combined to a dipole cross section (4.7). It is straight- 
forward to iterate this relation 



_ m 

C u 



m 



1 



dCn(C)<r(u-y) 



-ikj_-(u-y) 



(4.10) 



Here, the prefactor stems from the longitudinal posi- 
tion ordering £j < of subsequent interactions. The 



sum (4.6) over arbitrary m-fold rescattering contribu- 



tions can thus be represented by a simple exponential 
of the dipole cross section 



F (u,y) 



= P -* k i-( u -y) 



xexp^-i^ + d?n(0 a (u-y)^ . (4.11) 

2. Rescattering for yi < £, < yi 
The rescattering in the region yi < £ < yi is described 



by the path integral M of (4.2). To simplify this expres- 
sion, we start from its n = 1 first order contribution in 
opacity 



M^ 1 \r 1 ,v 2 ,r 3 \y h y l )=C A 



dqi± | / >,2 
[2 M<?i±)| 



(2tt) 2 



dCinfo) /iri(Si)dr 2 (fi)£iraKi) 
xMo(ri,r2,r 3 |y / ,^i)Mo(ri,r2,r3^i,2/ / ) 

x ^qii-(r 3 (6)-ri(?i)) _ A _ (4^) 

Here, Mo denotes the non-interacting part of M, ob- 
tained by replacing G — > Go in (4.2). The dipole cross 
section (4.7) can be identified with the qi^-integral. 
Equation (4.12) allows for an iteration to arbitrary or- 
der n. Writing the free Green's functions of Mq in a 
path-integral representation, using r = dr/d£, we find 



M(ri,r 2 ,r 3 |yi,yj) 



X>n Vr 2 Dr 3 



x exp 
x exp 



E f i - 
*T / <^("; 
z J m 

d£,n{£)a (r 3 - n 



1 - x) vl - f 2 



'HI 

1 f s! 

'24 



(4.13) 



3. Radiation cross section in terms of the colour dipole 

The medium-dependence of the radiation cross section 
(4.1) is contained in the functions F and M of ( 4.11 ) and 



(4.13). It is parametrized in terms of the colour dipole 
cross section a and the density of scattering centers 
Inser ting the expressions (4.11) and (4.13) for F and M 
into (4.1), a substantial amount of straightforward alge- 
bra allows for furt her s implifications. Two of the three 
path-integrals in ( 4.13| ) and some of the transverse in- 
tegrals in (4.1) can be done analytically. We give these 
technical details in app endix |5| The final result for the 
radiation cross section ( 3. 13| ) depends to leading O(x) on 
one path integral only 



lC(r(yi),yi;r(yi),yi\u) 



Vr exp 



r Vi ( I.) 1 



and takes the explicit form 
dV m ) a, 1 



d{\nx)dkj_ (2tt) 2 uj 2 



N c C F 2Rc 



(4.14) 



dyi / dyi 



-Avi\ -*\vi\ 



du e 



-ik 



d d 

x it ' Tr^y = Q,yi-, u ,yi\u) ■ 

oy ou 



(4.15) 



This is the main result of the present section. It contains 
all the characteristic features of the result of Zakharov [ 
H , and the extension of Zakharov's result to finite trans- 
verse momentum [ [l2[ [l3|. The dipole cross section cr(p) 
is the leading 0(x) contribution of the q-q-g cross sec- 
tion introduced by Nikolaev and Zakharov [ ||l[] and used 
in [ H |l3|] . The reason why this cross section depends 
on only one impact parameter p is evident from ( 3. 13 ) : 
in principle, rescattering contributions can compare the 
transverse position between the three combinations q — q, 
q — g or g — g in amplitude and complex conjugate am- 
plitude. In the leading 0(x) approximation, however, 
the transverse position of the quark remains unchanged 
due to rescattering. Moreover, the difference between the 
transverse paths of the gluon in amplitude and complex 
conjugate amplitude which is given by the function F in 
( 4 . 1 1 ) , is fixed by the measured transverse momentum 
and shows no time (i.e. £-) evolution. Accordingly, there 
is only one remaining impact parameter which measures 
the transverse separation of the gluon from the quark. 



The result (4.15) is closely related to the QED radia- 
tion cross section, see Eq. (5.13) below. The result (|4.15| ) 
differs from the expressions given by Zakharov [ ||, |l3| by 
(i) including th e e-r egularization and (ii) not subtract- 
ing from K. in ( 4.15 ) the zcroth order term /Co- As dis- 
cussed elsewhere in detail [ |l2| , (i) is indispensable for a 
quantitative analysis of the radiation cross section. The 
subtraction of /Co wa s included in Ref. [ |[ |l3| to cancel 
a potentially dangerous singularity of the path-integral 
/C. However, in the presence of the e-regularization, this 
singularity does not exist and the /Co-term vanishes. It 
can thus be included for technical convenience, but it is 
not needed. 
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Our derivation of ( 4.1 5| ) provides a new and direct 
proof of the equivalence of the BDMPS- and Zakharov- 
formalism. In fact, the non-abelian Furry approxima- 
tion which allowed us to write the high-energy limit of 
( [2.4|) is a compact shorthand of all the rescattering dia- 
grams included in the analysis of BDMPS [[!§]. Obtain- 
ing from this starting point the path-integral expression 
(4.15) proves the equivalence. It is an independent confir- 
mation of the corresponding statements in Ref. [ [l0| , and 
extends these arguments explicitly to the k^-diffcrential 
radiation cross section. 



V. PROBABILITY CONSERVATION TO FIXED 
ORDER IN OPACITY 

In this section, we introduce a new quantity a c \, ob- 
tained by factorizing out of the radiation cross section 
(4.15) a momentum independent absorption factor, 



dh 



d(lnx) dk± 



d 3 a r 



d(lnx) dk± 



(5.1) 



Here, Vtot denotes the total elastic cross section for the 
projectile scattering off one scattering center. In what 
follows, we argue that the iV-th order opacity term of a c \ 
reduces in the coherent and incoherent limiting cases to 
the classically expected results for bremsstrahlung asso- 
ciated to TV-fold scattering. In the incoherent limit, it 
becomes the incoherent sum of N Gunion-Bertsch radi- 
ation terms 



hm f^n 



L no— fixed 



V* NcGf ^\ 



N 




(5.2) 



In the coherent limit, it reduces to a single Gunion- 
Bertsch radiation term for the combined momentum 
transfer of N single scattering centers, 



1UO ^(AQ 

L^o d(\nx) dk 




a s AT n {L n ) 



N 



(5.3) 



In what follows, all limits L — > and L — > oo are taken 
at fixed opacity L«o = fixed without specifying this ex- 
plicitly. The N-th order contribution of (5.1) reads 



d 3 a^(N) 
d(ln x) dkj_ 



i -nm d 3 a c i(N-m) 
(-1) w v 



m=0 



d(ln x) dkj 



(5.4) 



where the weights w m are the m-th order terms of the 
absorption factor 



M'r, 



(5.5) 



These weights will be seen to appear as a consequence of 
probability conservation. In section |V A| , we discuss the 
above relations for QED. In section |VB| , we extend this 
discussion to QCD. 



A. Elastic scattering and radiation in QED 



In the present subsection, we take 
V(q±) 



«o(<Zx)| 2 
dqj_ 
(2nf 



ao(q±)\ 



(5.6) 
(5.7) 



to be the differential and total elastic scattering cross 
section for an electron scattering off a single electrically 
charged static potential. In general, we distinguish by a 
bar the QED quantities used in this subsection from their 
QCD counterparts discussed in the rest of this paper. 
The QED dipole cross section is defined by 



1-e 



(5.8) 



and the Fourier transform of a takes the form [ 



£(qx) = \ 



dp 



a{p) e lCil 



2tt) 2 

= -y(q ± ) + Vt ot ^ 2 )(q ± ). 



Furthermore, we denote by 



£(kx,q±) = 



2 2 



ki (kj 



-xqj_) 



(5.9) 



(5.10) 



the momentum dependence of a Bethe-Heitler radiation 
term associated with the incoming electron momentum 
kj^ and the scattering momentum transfer q^. Finally, 
we introduce two shorthands for the average over the 
transverse momenta q^j_ of the different scattering cen- 
ters. Depending on w hethe r this average involves the 
term cx Vtot ^Hqx) in ( 5 . 1 2| ) or not, we write 



= (-1)^11 

i=l 

dq.j± 
{2-k? 



dqix S(qix) 
ao(q*x)| 2 ■ 



(5.11) 
(5.12) 



The QED radiatio n cro ss section derived in [ |lj, [ll|] can 
be obtained from ( 4.15 ) by substituting the QED dipole 
cross section a of (5.S) for the QCD one, and by changing 
the coupling strength a s N c Cf — > OL exa , 
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d 3 a (in) 



d(lnx) dk± 



xe 



9 d rt 
oy ou 



due 



dyi / dyi 

J Vl 

-i df>(C)<x(u) 



0,yj;u,j/i|w) . 



(5.13) 



contrast to the QCD case ( 4.15| ), the term 
5 ^ n (^) ^( u ) ) m (5-13) ranges from z_ to y; [ 



In 

exp 

|j~2t since it is the rescattering of the incoming charged 
projectile (and not the rescattering of the outgoing ra- 
diated particle as in QCD) which determines the QED 
radiation spectrum. 

The following argument is based on the observation 
that the term oc Vtot 8^ 2 \<l±) in (5-9) can be factored out 



of (5.13) in form of the absorption factor in (5.1). The 
quantity a c \ depends only on V(q_i_). Moreover, the iV-th 
order term of the radiation cross section (5.13) depends 
only on the (Fourier transform) of the dipol e cro ss sec- 
tion (pj|), i.e., it is written with the averages ( 5.11 ). This 
suggests that a c \ tak es a simple form if written in terms 
of the averages ( pM2| ). To make more specific statements, 
we now inv estiga te the first few terms in the opacity ex- 
pansion of ( 5.13 ) which were calculated in Ref. [ 

The N = 1 term is the Bethe-Heitler radiation cross 
section times the opacity L uq of the medium, 



d(ln x) dkj_ 7r 



(Ln ) i?(kj_,qij_) 



-y(ijio) / i?(kj_,qij 

11 Jq.i 



(5.14) 



This confirms ( p.4[ ) up to first order in opacity. Since 
the Bethe-Heitler term vanishes for vanishing momentum 
transfer, the two medium averages (5.11) and (5.12) make 
no difference for this term. This is different for the higher 
order terms for which we discuss here the coherent and 
incoherent limits. 

For two scattering centers, the radiation cross section 
takes the coherent limit [ [l2j 

lim d3 ° (in) ( N = 2 ) 
l^o d(\n x) dk± 



a 



cm 



(Ln o y 



7T^ 



2 

(Ln ) 2 



#(k_L,qi_L +q 2 ±) 



i?(k ± ,qi ± + q 2 j_) 
qi J<i2 



J?(k_L,qi_ 



(5.15) 



Here, the first line is the result obtained in Ref. [ 
|l2|, the second line is obtained by inserting expression 
( |5.9|) and integrating out the ^-functions of the terms 
oc Vtot 3^ In the incoherent limit, we find in the 

same way 



d 3 a^(N = 2) 
l^oo d(lnx)dk± 
a em {Ln ) 2 



a cm (Ln ) 



[fl(kj.,qa ± ) 
R{k± + qi±,q 2 ±)] 
R(k± + qi±,q 2 



-Vtot / i?(k_L,qi_ 



(5.16) 



The results ( 5.15 ), ( 5.16) are the QED analogue of the 
expressions ( 3.47 ) and ( 3.48 ) obtained for t he N = 2 
QCD radiation cross section. This confirms (5.4) up to 
second order in opacity: 

d z a^ n \N = 2) rf 3 cr cl (7V = 2) d 3 a cl {N=l) 



d(lnx) e?kj 



d(]nx) dkj 



Wl 



d(ln x) dkj 



(5.17) 



Here, the weight w\ = (Lno Vtot) contains only informa- 
tion on the mean free path of the projectile. It is the 
probability of having one additional interaction with the 
medium. The physical reason for the negative terms in 
( |5.15 ) and ( 5.16Q is now obvious: 
An expansion of (5.13) up to second order sums over 



the two possibilities that the photon was produced by an 
interaction of the electron with either one or two scatter- 
ing centers. The probability of this second interaction is 
a factor w\ smaller than the probability that only one in- 
teraction occurs. Hence, to go consistently from the first 
( N = 1) to the second (N — 2) order approximation of 
( p.l3|) requires two steps: (i) the radiation contribution 
related to two rescatterings of finite momentum transfers 
has to be added [this contribution is denoted by a c \] . (ii) 
the probability for the N — 1-contribution has to be re- 
duced by a factor (1 — w±) [in this way, the probabilities 
for the N = 1 and N = 2-cont ribut ions add up to unity] . 
The second, negative term in (5.17) implements step (ii). 
It is this second step which leads to the negative term 
in the coherent limit (5.15) and which changes the sign 
of the second term in the incoherent limit (5.16). To 
sum up: the opacity expansion of (4.15) correctly read- 
justs the probability of the cr c \(N — 1) contribution if 
including the higher order = 2 in the approximation 
of ( gl| ). 

Ref. [ |l2| also provides the third order contribution to 
the QED radiation cross section. The limiting cases read 

lim d3 ° (in) ( N = 3 ) 
l^o d(lna;)(ikj_ 



0,. 



(Ln f 



3! 

-3 V tot 



i?(k|_,qi_L + q 2 ± +q3J 



qi J<12 Jq3 



q 2 j 
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R(k±,cn ± ) 



(5.18) 



r rf35 ' (m) (^ = 3) 



Q » AT" ^ n 2 ( Ln 0) 

— l\ c Of O a 



L no— fixed 

2 J R(k L +gu_,qr 2 x) + i?(k ±! 9l± )] 



91 J 92 J 93 



3! 



(5.19) 



The relative signs and weights of different Bethc-Hcitlcr 
terms confirm (5.1) to third order 



d 3p(in) (JV" = 3) _ d 3 a cl (N = 3) d 3 a cl (A^ = 2) 

ci(ln a;) dk± c?(ln x) dk± d(ln a;) dk±_ 

d 3 a cl (N = 1) 



+w 2 - 



d(ln x) dk.± 



(5.20) 



We note that summed over arbitr ary orders of opac- 
ity, the factor X^o( — ^ w i m (^i) am ounts to a 
momentum-independent normalization. The angular dis- 
tribution of gluon ra diat ion, calculated directly from the 
all-order expression (5.1) is not modified by this normal- 
ization factor. In an expansion in opacity, however, the 
weights Wi do not factorize from the momentum depen- 
dence. The correct angular distribution for radiation off 
a target of N scattering centers is given by the sum over 
the opacity contributions up to N-th order. In this way, 
the results of the BDMPS/Zakharov- formalism can be 
compared directly to results for a finite number of scat- 
tering centers. 



B. Elastic scattering and medium-induced radiation 
in QCD 



The above arguments confirm the expressions fl5.l| )- 
( |5.5D for QED, and they apply to the QCD radiation 
spectrum, too. This follows from the fact that the QED 
and QCD radiation spectrum are related by s imple sub- 
stitutions which do not affect the form of ( 4.15 ). In QED, 
however, the only total elastic scattering cross section 
Vt t associated to the projectile is that of the electron. 
In contrast, both quark and gluon are charged in QCD 
and have a finite elastic scatteri ng c r oss section. What 
is additionally needed to justify (5T)-(5 : 5) in QCD is an 
argument why the total elastic scattering cross section 
Vt t in (5.1)-( |5~5| ) should be that of the gluon, 



C A 



dqi 

(2tt) 



M<Z_l) 



(5.21) 



rather than that of the qua rk or a combination of both. 

The mot ivatio n for ( 5.21 ) comes from the observation 
of section III A that the i-th scattering center always 



contributes with a factor Ca to the radiation cross sec- 
tion, i.e., the (?-<?-projectile scatters effectively with the 
coupling strength of a gluon. This is a consequence of 
the leading 0(x) approximation in which, as observed by 
BDMPS, only the rescattering of the gluon matters. Be- 
yond the leading 0(x) approximation, more complicated 
combinations of elastic quark and gluon cross sections 
may be needed, but this is the regime in which colour 
triviality breaks down and our formalism does not apply. 

I n an alogy to th e QED case, we can argue on the basis 
of ( |5.l| ) and ( 5.21 ) that the A^-th order opacity term of 
the QCD radiation cross section (1.15) is a convolution 
of the A^-fold scattering contribution and a readjustment 
of the probabilities that less than A?" rescatterings occur. 
This explains especially the structure of the N — 2 QCD 
contributions calculated in section III Ej . More details 
about the opacity expansion for QCD are given in the 
next section and in Appendix |^. 



VI. IN-MEDIUM PRODUCTION: HARD VERSUS 
MEDIUM-INDUCED RADIATION 



The initial state boundary condition in the above cal- 
culations is tha t of a quark approaching an asymptotic 
plane wave (2J3) at far backward positions. In rclativistic 
heavy ion collisions, however, the hard (off-shell) parton 
is produced inside the collision region and its propaga- 
tion inside the medium has to be followed from a fixed 
finite time onwards. In what follows, we discuss the 
corresponding gluon radiation cross section a^ nas ^ for a 
nascent quark, in contrast to the case for a free incoming 
quark which was discussed in sections [n|-0. 

The radiation cross section a^ nas ^ for nascent quarks 
is obtained by making two simple modifications to our 
derivation in sections ^-|v|: 

(i) At the beginning of our derivation in section [n|, one 
has to replace the incoming plane wave -F(x,pi) by a 
production amplitude J{pi) which specifies the probabil- 
ity that at some point (z , 0^, z ) in the collision region 
a hard parton of momentum pi is produced. We follow 
Gyulassy, Levai and Vitev [ [l7j] in factorizing the pro- 
duction probability |J(pi)| 2 out of the radiation cross 
section. This is an assumption on the weak momentum 
dependence of this source which is also implicitly present 
in other discussions of radiative energy loss from partons 
produced in medium [ ^j. Since our starting point (2. IS) 
is written for pi^ = 0, i.e., -F(x,pi) = 1, this amounts to 



multiplying ( 2. IS ) by a factor | J(pi)| . To simplify nota- 



tion, we drop this factor in what follows: in comparing to 
experiment, the radiation cross section cr( nas ) given be- 
low has to be weighted with the production probability 
of the hard parton. 

(ii) The parton produced at (zo, 0^, zq) propagates in the 
positive z-direction. Accordingly, the gluon emission ver- 
tices in the amplitude and complex conjugate amplitude 
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are constrained to longitudinal positions yi, yi > Zq. This 
changes the boundaries of the corresponding integrals in 
( pl8| ) and survives all interm ediat e steps. 

With these two changes in ( 2.1 8| ) , and choosing zq = 0, 
the medium-induced gluon radiation cross section off a 
nascent quark becomes 



d(hix) dk_ 
x2Rc 



1 



N g Cf 



\vi\ 



(2tt) 2 lo 2 

Z+ ,2+ 

dyi / dy t e-^' 

-i P+ «n(Oa(») 



due" lk - 

" ^"^(y = 0,yj;u,y/ a;) 
ay ou 



(6.1) 



This differs from the corresponding result of Zakharov [ 
|l3f by the regularization prescription onl y. 

In what follows, we shall analyze (6T) in the opacity 
expansion. To this end, we expand the path-integral /C 

z 

Z 

z' 

+ J d£in(£i) J d&n(&) j dpi dp 2 /C (r, yi; p u fi) 

xifT(p 1 )/C(pi,^i;p 2 ,C2) ^ cr (/°2)^o(P2,6;r,y/) ■ (6-2) 

Here, we have suppressed the explicit w-dependence in 
/C. The corresponding free Green's function /Co reads 



fco(ir,yi;r,yi) 



2ni(yi - y t ) 



exp 



iu> (r — r)' 
2(Vi-Vl) 



(6.3) 



The technical steps involved in the opacity expansion are 
described in detail in Ref. [ [l2) . We present our result for 
the N-th. order opacity contribution to (|6 . 1[) in the form 



d 3 a^(N) _ a s 1 
d{\nx)dk ± (2tt) 2 uj 2 C F j- 



N 



N 



A*) 

-N.J+l 



1=1 




(6.4) 



where the integration is defined by fl5.ll ) and the fac- 
tors Zftj are tabulated in Appendix |j. The superscript 
(* = in or * = nas) indicates that the result differs for 
the radiation cross sections ( 4.15 ) and (3.1) by the spe- 



cific form of the factor Z K N '■ only. We have checked the 



form of (6.4) for N < 3 by explicit calculation, but we 
expect that it holds for arbitrary N. Especially, it is ob- 
vious that the N-th. order term contains exactly (N + 1) 



different terms: the j-th term in ( |6.4| ) corresponds to 
the (N — j)-th order expansion in /C and the j-th or- 
der expansion of the exponen tial exp [— | J d£ n(£) cr(u)] . 
Moreover, each integrand in (6.4) has to de pend on all N 
transverse momenta, since the integration ( |5.11 ) leads to 
a vanishing contribution otherwise. 



The N-th order opacity expansion of (6.1 ) involves only 
(N + 1) different terms. This is a significant simplifica- 
tion in comparison to brute force calculations [ [l7| which 
require already for TV = 2 the calculation of 7 2 direct plus 
2 x 96 contact terms [ ^3| . In what follows, we discuss the 
physics contained in this opacity expa nsio n by studying 
the coherent and incoherent limits of (6.4). 



A. Hard radiation and its N = 1 modification 



For the gluon radiation spectrum ( 4.15| ) of a free in- 
coming quark, the zeroth order in opacity vanishes. A 
free incoming particle does not radiate without interac- 
tion. A nascent parton, however, does radiate without 
further interaction. To zeroth order in Luq, we find 



d 3 a (na8)( N = j as 



N C C F 



1 



d(\nx)dk ± 7T 2 ' " ~* k\ ' 
This is the characteristic gluon radiation spectrum 



H(k A 



1 

ki 



(6.5) 



(6.6) 



associated to the production of a hard parton. In energy 
loss calculations which integrate over u>, the above radia- 
tion cross section (3.5) as well as its medium-dependent 
form should be multiplied with the quark-gluon splitting 
function i _Hlz_) w hich equals one only in the O(x) ap- 



proximation [jl4|]. Equation (6.5) is the QCD analogue 
of the QED radiation accompanying /3-decay. To deter- 
mine the energy loss of hard partons in relativistic heavy 
ion collisions, we are interested in the medium-dependent 
changes of the radiation cross section (6.5). To first order 
in opacity, we find from Appendix [|] 

d 3 a( nas \N = 1) _ 



dQjax) dkj 



a. N C C F 

r2 



2iJ 2 



QQi 



xtiq ■ 



LQi - sm(LQi) 



(6.7) 



To understand the physical meaning of this expression 
it helps to rewrite the coherent and incoherent limiting 
cases in terms of the hard and Gunion-Bertsch term: 



lim 



d 3 a (nas)( N = y 

L—>h d(lnx)dk± 

d 3 (j( Ms )(jv = i; 



lim 



d(lnx) (ikj 



= 0, 



= -L N C C F (no L) 



[H(k ± + qix) + iZ(k X) qix)] 



(6.8) 



(6.9) 
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In the incoherent L —> oo limit, the radiation cross sec- 



tion (6.1) expanded up to first order in opacity, takes the 
form 



lira 



N=l 



dV nos ) (m) 
d(ln x) dk^ 



= -4 N c C P 



m=0 



wi)i?(k ± ) +n i / H(k ± +m_ 

J q! 



(1 



n L J R(k±,qx±) 

qi 



(6.10) 



The three terms on the r.h.s. of ( 6.1C ) have a sim- 
ple physical meaning: the first is the hard, medium- 
independent radiation fl6.6| ) reduced by the probability 
wi = uq L Vtot that one interaction of the projectile oc- 
curs in the medium. The second term describes the 
hard radiation component which rescatters once in the 
medium. The third term is the medium-induced Gunion- 
Bertsch contribution associated with the rescattering. In 
general, the L-dependence of ([T?]) reflects the interfer- 
ence pattern between the different contributions and can 
be quite complicated. For the case N = 1, we have seen 
that this interference pattern interpolates between simple 
and physically intuitive limiting cases. 



B. Coherent and Incoherent Limits for N = 2 and 

N = 3 

The general L-dependent expressions for the second 
and third orders in the opacity expansion of ( |6 . l| ) can eas- 
ily be written down by using ( |6.4| ) and the corresponding 
results tabulated in Appendix [El In the present section, 
we focus entirely on the coherent and incoherent limiting 
cases. For N = 2, they read 

d 3 {nas)( N = 2 ) 

lim — 771 L = > 6 - U 

l^o d(\nx)dk± 



.. riV" as >(/V = 2) a s (n a L 
um — r-^ = — o JVc Of 



d(lnx)dk± it 2 "" J 2 
H(k ± + q 1± + qax) + R(k x + qix.qax)] ■ (6.12) 



Adding this second order contribution to ( 6.1 0| ) , a simple 
systematic starts to emerge: 



N=2 



lim >^ 



d 3 ff'"° s »(ra) 



L—nya ' d(lnx)dki 

m=0 y ' 



= ^N C C F 



[{\-w 1 +w 2 )H(k 1 _) 

+(l- Wl )n L I H{k ± 
Jqi 



lim 



d 3 a cl (2) 
d(ln x) dkj 



H(k A 



+ (l-wi) 



■qi-U 
+- qi± + q2_ 

rfVcl(l) 

d(\nx) dk± 



(6.13) 



Both the m-fold rescattered hard term and the Gunion 
Bertsch contribution for m-fold rescattering are weighted 
with the expansion of the absorption factor up to order 
(N — m). For N = 3, we find from the results in Ap- 
pendix Fj 



d 3 a (nas)r N = 3) 

km — — — — = , 

l^o d(mx)dk± 

r rfV"^)(jV = 3) a s (n Lf 
lim — — = ~ Nc C f 



(6.14) 



i^^oo d(lnx)dk± 7T- 
x [H(k x + qi± + q2X + Q3x) 

J En 



3! 



-i?(k ± + q i± + q 2 _L, q 3 _L)] 



(6.15) 



Adding this contribution to ( 6.13 ) confirms the indicated 
systematics. Rather than spelling out the result, we ex- 
trapolate in the following section our analytical findings 
to arbitrary TV. 



C. Extrapolation of limiting cases to arbitrary N 

The results of our explicit calculations up to third or- 
der, combined with the physical arguments given above, 
suggest results valid to arbitrary orders N. For the case 
of a fr ee i nco min g quark, these results are given in equa- 
tions (|5.2| ) - ( 5A ) . For the case of an in- medium produced 
quark, we write them in terms of the shorthand 



H (m \k ± ) 




(6.16) 



This describes the hard radiation component which 
rescatters on m well-separated scattering centers. For 
the incoherent limit of ( |6.1[ ), expanded up to TV-th order, 
we obtain by extrapolating the systematics observed up 
toiV = 3: 



N rfZ^inas)/^ N 



N- 



L^oo ^ d(hix)dk\ 

m=0 y ' m=0 



N 



N-m 



WN~ 



lim 



w N - m H(' n Hk ± ) 
d 3 a c \(m) 



?TI = 



L— »oo <i(ln a;) dk± 
In the coherent limit, we infer from ( 6.11 ) and ( |6.14| ) 



(6.17) 



N 



d 3 a (nas)r ) X 

m— 



In this coherent limit, there is no medium modification to 



the hard radiation (6.5), irrespective of the opacity of the 
target: locating the medium at (and only at) the creation 
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point of the parton, the projectile looses the possibility 
of interacting with the medium. 

What matters physically is, how medium-effects switch 
on with increasing thickness L of the medium. This in- 
dependence will interpolate smoothly between the totally 
coherent and totally incoherent limits discussed here. 

VII. CONCLUSION 



incoherent limit confirms the classical picture of an inco- 
herent sum of rescattering contributions of the hard ra- 
diation and Gunion-Bcrtsch radiation terms. The depen- 
dence of gluon radiation on the thickness L of the medium 
is expected to interpolate smoothly between these ana- 
lytically accessible limiting cases. The detailed study of 
this L-dependence awaits further work, not only in the 
opacity expansion discussed here, but also in the dipole 
approximation discussed previously [ [l2| . 



Colour triviality is an important tool in the study of 
the nuclear dependence of hard QCD processes. For 
colour trivial observables, the A^-fold gluon exchange 
with the nuclear environment reduces to the iV-th power 
of a colour Casimir operator. This renders the calcula- 
tional problem abelian. 

In the present work, we have given a simple diagram- 
matic proof of the colour triviality for the iV-th order 
opacity contribution to the kj_-differential gluon radia- 
tion cross section off a hard quark. The diagrammatic 
identities developped to this end are a direct consequence 
of the non-abelian Furry approximation for hard partons 
in a nuclear environment. Summing up the contributions 
to the radiation cross sections to all orders in opacity, we 
have derived Zakharov's path-integral formalism and we 
have given a new proof of its equivalence to the BDMPS- 
formalism. 

For practical applications, it is important that our ap- 
proach gives an easy access to the low opacity expansion 
of the radiation cross sections ( 4.15 ) and (S.l). The iV-th 
order term involves only [N + 1) terms which are rela- 
tively easy to calculate. In contrast to previous brute 
force calculations [ [TtJ which have to deal with an ex- 
ponentially increasing number of terms, our approach is 
thus well-suited for the calculation of a realistic number 
of scattering centers. This is illustrated by the simple 
form of our explicit results up to third order in opacity. 

As discussed in section the information contained 
in the iV-th order opacity term to the gluon radiation 
cross section is a convolution of (i) the radiation asso- 
ciated to the rescattering off N scattering centers and 
(ii) a readjustment of the probabilities that rescattering 
occurs with m < N scattering centers only. We have 
shown how to disentangle this information by factoriz- 



ing in (5.1) the radiation cross section into two terms. 
In this way, at least the coherent and incoherent limits 
of the involved L-dependent interference pattern became 
accessible analytically to arbitrary order. 

Most important for applications to relativistic heavy 
ion collisions at RHIC and LHC, we have extended these 
results in section ^ to the gluon radiation off quarks pro- 
duced in the medium. The corresponding radiation cross 
section (6.1) gives access to the interference pattern be- 
tween the hard radiation (associated to the production 
of the hard parton) and the medium-induced radiation. 
Its coherent limit shows no medium dependence. The 
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APPENDIX A: NON-ABELIAN FURRY 
APPROXIMATION FOR RESCATTERING 
GLUONS 

In this appendix, we simplify the gluon rescattering 
contribution (2.7). We consider on-shell gluons = 
(u,kj_, ki) with a physical transverse polarization 



t(k) = (e ,e_L, -e ) . 



eo 



w + ki 



(Al) 



To leading order in energy, the contraction of this polar- 
ization vector with the momentum dependent part of the 



three-gluon vertices occuring in (2.7) takes the form 



V^iO/^Cfci-ijfei - h-1, h) e M, (fc) 
= -2we^- 1 (fc) + 0(cj°). 



(A2) 



This allows to rewrite (2.7) 



I^W(y,k)=e^(k)V (n/yfj 



-d 3 x- 

\3 a Xl 



x[-i4 9) W] 



d 3 k 

>) 

-i (-2w) 
kf + i e 



x gikj'(xj— Xj_i) \ g-ik-xt 



(A3) 



where we have rearranged the phases, using x = y. The 
longitudinal momentum integrals can be done by contour 
integration 



dk 



u 



(2tt) kf+ie 
1 



2 ki. 



Q(x itl -x i . ltl )e ik ^^'- x *-^ . (A4) 



On the r.h.s. of this equation, kij is determined by the 

,2 

pole value to order 0(l/w), fey = lu — -J^. Keeping in 
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(a.3) the norm to leading order in 1/u but the phase to 
next to leading order, we write 

J Ml(L) (y,fc) =^e- iun V (J] / d^eix^-Xi^) 



xe 



(A5) 



where the transverse momentum integrals are identified 
with the free Green's functions 



f trk, 
GoCxi-ijXilw) = J 



pikt.i'Cx^x— 



xe 2 



(A6) 



To rewrite the path-ordered product in (A5), we intro- 
duce the shorthand 

G {L) (y±,yi;x±,xi\u) = 

x ^JJ[-«A)(xi)] Go(xi,x,^i,i;xj +li x,^i+i,ib) J i ( A7 ) 
where x^ + i = x. This allows us to write 

x J dx ± G( L \y;x\iu)F(x,k). (A8) 



APPENDIX B: SIMPLIFYING THE 
INTERACTION VERTEX 

Here, we dis cuss details of how to simplify the inter- 
action vertex (2.16) in the leading 0(x) approximation. 
The y-derivatives in the differential operators D\ and D 2 
are conjugate to the transverse momenta pj\ and p 2 \ of 
the quark entering and leaving the radiation vertex, see 
Fig. [l] for notation. As a consequence, the interaction 
vertex ( |2.15D takes a particularly simple form in a frame 
in which pj^ = p^i . Depending on the spin of the ingoing 
(X q = ±o) and outgoing (X q / = ±\) quark and the gluon 
helicity (X g — ±1), it reads in this frame (see Refs. [ 
pd| , for a derivation) 



Ty(X q — Xg',Xg) 



i X g [X g (2- x) + 2X q x] 

— -iX — 

dyi 9 dy 2 



r r ( Xn — —X„' , A„ 



= 2m q xX g 5\ g ,2A, , 



(Bl) 
(B2) 



where the differential operators in (Bl) act on the two dif- 
ferent transverse components y = (3/1,2/2) of the Green's 
function for the outgoing quark. The important property 
of this representation of the spinor structure is that the 
spin- and helicity-averaged combination f y f ^, takes the 
simple form 



ft - A_ JL 



■9sf-. 



(B3) 



where the prefactors for the spin-flip g s f and non-flip g n f 
term read 



g nf = [4 - 4x + 2x 2 } . 



9s f = 2m 2 q x 2 



(B4) 



For the case of photon emission, one can choose to rede- 
fine the longitudinal z-axis to be parallel to the emitted 
photon k [0,|l2|. This ensures pj\ = p 21 and allows to 
eliminate all spinor dependence from the radiation proba- 
bility with the help of ( p33| ) and y* = y. For the gluon ra- 
diation amplitude, this trick cannot be used. The reason 
is that due to the rescattering of the gluon, the direction 
k of the final state gluon is not aligned to the direction 
ki with which the gluon leaves the radiation vertex. To 
use the simple expression (|B3| ) nevertheless, we introduce 
the derivative —id/dy g , where the subscript denotes that 
this partial derivative is acting on the y-dependence of 
the gluon Green's function G( s ). —id/dy g is the conju- 
gate of the direction ki j_ . In the frame parallel to ki , the 
conjugate of the momentum p^\ thus takes the form 



d_ 
dy 



d 
dy* 



JL 
dy 



l-x d 
x dy g 



(B5) 



With this replacement, the vertex functions in (Bl) and 
(B3) can be used for the discussion of the gluon radia- 
tion spectrum in a frame in which the incoming quark 
momentum pi defines the longitudinal direction. 

In the leading 0(x) approximation, only the partial 
derivative w.r.t. y g survives in (B5), the spin-flip contri- 
bution can be neglected, and the spin- and helicity av- 
eraged combination (B2) takes in momentum space the 
simple form ( 2.25| ). 



APPENDIX C: CANCELATION OF A CLASS OF 
CONTACT TERMS 



To arrive at the expression (3.13) for the gluon radia- 
tion cross section in the presence of N scattering centers, 
we have included arbitrary combinations of real and con- 
tact terms to fixed order in opacity. However, contribu- 
tions to the radiation cross section have to involve at least 
one (real or contact) interaction in both the amplitude 
and complex amplitude. The reason for that is that due 
to energy-momentum conservation, an amplitude which 
is touched neither by real nor by contact terms cannot 
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contribute to a final state. According to this counting, 
diagrams which have N contact terms in one amplitude 
should not be included in the calculation of the O(N) 
contribution. Here, we show that their inclusion in the 



derivation of ( 3.13 ) is allowed, since this set of diagrams 
adds up to zero at arbitrary order in opacity. 

For our proof, we consider the set of O(N) diagrams 
with N contact terms in the amplitude, and m of these 
contact terms occuring at longitudinal positions > yi. 
These contributions have the form 



N-m m (CI) 

where each thin vertical line denotes a contact term link- 
ing either twice to the quark line or twice to the gluon 
line or once to the quark and once to the gluon line. In 
what follows, we consider the subset of these diagrams 
which contributes to an arbitrary but fixed interaction 
vertex rVj). The longitudinal phase factor <!> associated 
with this set of diagrams takes the form 



®(yi,yi) = <p{yi) 



-iQyi 



(C2) 



The reason for this factorization of the ^-dependence is 
that the longitudinal positions £», i € [1, AT] , of the N 
scattering centers (which are integration variables in the 
calculation of $) do not have yi as integration boundary. 
As a consequence, the only possible yj-dependence of $ 
comes from the phase in the conjugate interaction free 
amplitude where the gluon has transverse energy Q. 

For each class of diagrams (CI), we find the corre- 
sponding class of diagrams which has N contact terms 
on the opposite side of the cut 



N-m (C3) 

For the class of diagrams to fixed rVj), the phase corre- 
sponding to (C3) reads 



®*{vuvi) = <P*(yi)e 



(C4) 



It is related to ( p2| ) by complex conjugation and inter- 
change yi <->• yi. The combined Z- factor for the sum of 



both phases (C2) and (04) reads according to (3.19) 

■^comb 



2Re f + dyi f + dyi e - e l«l- e l*«l 
Jz- Jyi 

x ($(yi,yi) + $(yi,yi)) ■ (C5) 



The second term in ( jC5| ) can be rewritten in the form 



dyi / dyi®{yi - yi) HvuVi) 



z + 

dyi I dyi 6 {y t - yi) $(yi,yi) . 



(C6) 



This allows to write for ( CJ ) 

dyi J dyi 



(C7) 



The yi-integration of this expression can be done explic- 
itly and results in -Z C omb = 0. This proves our claim that 
the sum of all O(N) contributions which have N contact 
terms on one side of the cut, vanishes. 



APPENDIX D: DERIVATION OF THE 



RADIATION CROSS SECTION (4.15) 



The radiation cross section fl4.15|) is obtained by in- 

"1 



sorting the path-integral M of ( 4T3h and the function F 
of (111) into ( 4.11 ). To simplify the path- integral M, we 
change coordinates 



r a (£)=xr 3 (0 + (l-z)r 2 (0, 
r 6 (0=r 3 (0-r 2 (0, 



(Dl) 
(D2) 



and we introduce relative and average pair coordinates 



0(O = r«(O-'i(O. 
p(0 = r»(0+ri(0- 



(D3) 
(D4) 



In terms of these coordinates, the path-integral M of 
CTl reads 



M(ri,r 2 ,r 3 |!/|,i«) = J VpVpVr b 



x exp 
x exp 



E f m 

*y J d( (p-p+ {l-x)xv h 



1 



<tt,n{£)a (p + (1 - x) r b ) 



(D5) 



The path-integral over Vp can be done analytically (see 
Appendix B of Ref . [ |l| for more details on this step) 



M(ri,r 2 ,r 3 |y z ,y z ) 



E 



s 2ni(yi - yt) 

n ,J E , [ 2 y-(v+s)+v 2 -s 2 l 



x / T>r b exp 



1 



/// 



dl; [i 



Hi 



.E(l -x)x. 



^ n (0 a (Ps + (1 - x) r b ) 



(D6) 



where 
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v = x(u - y) - xr 3 (yi) - (1 - x)y , 
and p s denotes the straight line path 



PsiO =x(u-y) - (v + s)- 



yi - yi 



(D7) 



(D8) 



Here, we have kept the dependence on r 3 (?/;) explicitly, 
since we need the corresponding derivative in (4.1). We 
now shift the intergration variable u — > u + y in (4.1), 
replacing the partial derivative w.r.t. y by the partial 
derivative w.r.t. u and doing a partial integration: 



I' 



A± 



dy dy ds du F(u) 



du 



8 



dr 3 (yi) 



M^+y\v 1 ,v 2 ,v 3 \y l ,y l ). 



(D9) 



Here, the superscript M ( u ^ u +y) on j\ f in dicates that the 
u variable is shifted. Starting from (D6), the y- and s- 
integration is now trivial: 



J dyd8M( u -* u +y)(ri,r 2) r3| W ,w) 



Vv b exp 
1 



-n(£)a (iu + (1 - x) r b ) 



(D10) 



The approximations used to derive (4.1) are valid to lead- 
ing order 0(x). To this order, the path-int egral in (D1C) 
coincides with the p ath-in tegral K, given in ( 4.14 ) . To ob- 
tain our final result ( 4.15 ) for the radiation cross section, 
we note that the boundary condition on this path-integral 
is r b (yi) = r 3 (yi) - y = and r b (yi) = u. We can thus 
replace in (D9) -g— -, — r — > -g- - ( — r. The integrand of (DS) 



9r 3 (j/j) dr b (yi) 

is then y-independent, and the y-integral cancels against 
the total transverse area A i . This leads to the result 



(4.15) 



APPENDIX E: GLUON RADIATION UP TO 
THIRD ORDER IN OPACITY 

Here, we give explicit results for the low opacity ex- 



pansion of the gluon radiation cross sections (4.15) and 
(3.1) up to third order. We have assumed a medium 
of homogeneous density no and finite thickness L. By 
explicit calculation, we have found for N < 3 that the 
iV-th order terms are of the form (6^). We denote by 
Qj the transverse energy of the gluon after j momentum 
transfers from the medi um, see ( 3.3£ ) for definition. The 
Z-i actors entering (OA) are: 



Results for 



dV'">(iV) 

d(\i\ x) dkj^ 



A. For N = 1: 



7 (in) 



7 (in) _ 
^1,2 — _ 



2Q\ ' 
Ln 



B. For N 



Z 2,l — 71 



7 {in) 
J 2,2 



^2,3 — n ( 



°4Q1' 
2 2 - 2cos(LQi) - L 2 Q\ 
2QI Q 2 

-1 + cos (LQx) 

QQl Q2 



C. For N = 3: 



(in) _ 

3.1 — 

(in) _ 

3.2 — 

7 (in) _ 
^3,3 — 



7 (in) _ 
^3,4 — 



12 Ql ' 

3 6sin(£Q 2 ) -6LQ 2 + L 3 Qj 
~ U ° 6QtQ 3 
3 Q\ [sin(£Q 2 ) - LQ 2 ] 

l ° Q\{Qi-Q2)QlQ3 

3 Qj [sin (LQi)-LQi] 

~ n ° g?(Qi-Q 2 )g|Q 3 ' 

3 Ql [LQ 2 - sin (£Q 2 )] 

l °QQj (Qi-Q 2 )QlQs 

3 Q|[LQi - sin (LQi)] 
Q Qi (Oi-Qa)Ql Qa' 



Results for ^ - . 



A. For N = 1: 



7 (nas) 



n 



7 (nas) 
^1,2 



LQi - sin (LQi) 
LQi - sin (£Qi) 



-n ■ 



B. For N = 2: 



Z 2,l 



7 (nas) _ 2 
'2.2 — ''0 ' 



2 2cos(LQ2) - 2 + L 2 Ql 
*° 2Q| ' 

Q? [2 -2 cos (LQ 2 ) -L 2 Q|] 



^(nas) 
J 2,3 



2Q\ (Qi - Q 2 )Q\ 
3 Ql [2-2cos(LQ 1 )-L 2 Qf] 
l ° 2Q\{Q 1 -Q 2 )Ql 
Q\ [-l + cos(LQ 2 )] 



QQl (Qx-Q 2 )Q 2 2 

3 g| [-l+cos(LQi)] 
'° QQl{Qi-Q2)Q 2 2 



(El) 
(E2) 

(E3) 
(E4) 
(E5) 

(E6) 
(E7) 

(E8) 

(E9) 



(E10) 
(Ell) 

(E12) 
(E13) 
(E14) 



C. For N = 3: 
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y{nas) 
J 3,l 



J 3A 



3 6sin(LQ 3 ) - 6LQ3 + L 3 Q 3 3 



6QI 



Anas) _ 3 
J 3,2 — "0 



J 3,3 



Q| [6 sin (Lga) - 6LQ2 + £ 3 Q%] 
6Q!(Q 2 -Q 3 )Qt 

3 Qj [6 sin (LQ 3 ) - 6LQ 3 + L 3 Q 3 3 ] 

'° 6g|(g 2 -Q 3 )Q! 

n 3 , sin (LQi) 



Ql (Qi - Q2) (Qi - Q 3 



ng sin (LQ 2 ) 



Ql (Qa - Qi) (O2 - Qs) 
n| sin (LQ 3 ) 

'Qi(Q 3 -Qi)(Q3-Qa) 
r Q 2 i93 + QiQ2 + QiQ3 



Qi Ql Ql 

ng sin(LQi) 



QQ?(Qi-Q 2 )(Qi-Q 3 ) 



rig sin(L<5 2 ) 



Qi (Q 2 -Qi)(Q 2 -Q 3 )Q 



sin (LQ 3 ) 



Q3 (Q 3 -Qi)(Q 3 -Q 2 )Q 



QQ1Q2 Q 3 ' 



(E15) 



(E16) 



(E17) 



(E18) 
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